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We argue that complex systems, defined as non-Poisson renewal process, with complexity index 
fi, exchange information through complexity matching. We illustrate this property with detailed 
theoretical and numerical calculations describing a system with complexity index fis perturbed by a 
signal with complexity index fip. We focus our attention on the case 1.5 < fis < 2 and 1 < jip < 2. 
We show that for fis > fip, the system S reproduces the perturbation, and the response intensity 
increases with increasing fip. The maximum intensity is realized by the matching condition fip — fis- 
For \ip > lis the response intensity dies out as l/t MP-Ms . 

PACS numbers: 05.40.Fb, 05.60.Cd, 02.50.Ey 



Aperiodic stochastic resonance lj and phase synchro- 
nization are recently discovered phenomena through 
which physicists try to establish a control on complex 
processes, notably those of neuro physiological interest. 
We study an intermittent behavior, reminiscent of the 
noise-free stochastic resonance process through chaotic 
maps -'i. |4|. However, we go beyond the choice of a driv- 
ing force either of Poisson kind or totally random. As 
to the synchronization, we establish it using the concept 
of Complexity Matching (CM). Our departure point is 
the recent observation [5J, |6j that a non-Poisson renewal 
system with power index fis < 2 does not respond to a 
harmonic perturbation. We devote this letter to proving 
that a non-Poisson renewal system with index fis < 2, 
on the contrary, is so sensitive to non-Poisson renewal 
perturbations with power index fip < 2 as to inherit in 
the long-time limit their power indexes. Furthermore, 
we shall show that the response intensity is maximum at 
the matching condition fis — fip, and that, when depart- 
ing from it, the response intensity either will be weaker 
{lis > /ip) or will die out (/is < fip). 

The CM phenomenon is a property of Non-Poisson Re- 
newal (NPR) processes, here considered as a paradig- 
matic example of complexity. The Blinking Quantum 
Dots (BQD) are NPR systems Q: The jump from 
the "on" ("off') state to the "off' ("on") resets to zero 
the system's memory, thereby ensuring the renewal con- 
dition. Furthermore, the time distance between two con- 
secutive events is described by a histogram with the form 
of an inverse power law with power index fis < 2. The 
free dynamics of the process are described by the dichoto- 
mous variable £,s{t), with the values £s = 1 and £s = — 1 
indicating that the system is in the "on" (or |1)) and 
"off' (or |2))state, respectively. As pointed out by Bel 
and Barkai the condition /is < 2 of the NPR sys- 



tems generates ergodicity breakdown, a condition shared 
by the phenomenological models of glassy dynamics |l0| , 
laser cooling [ill and models of atomic transport in op- 
tical lattices |l2j. 

To ease the theoretical treatment, we make the simpli- 
fying assumption that the "on" waiting time distribution 
is identical to the "off' waiting time distribution. When 
the system moves from an initial out of equilibrium con- 
dition, the survival probability ^s(t) is given by 



*s(t) 



t + T s 



(1) 



This experimentally observable relaxation is traced back 
|l3| to the waiting time distribution 



*Ps(t) = (us ~ 1) 



rpfJ.S-1 

(r + TsY 



(2) 



This distribution refers to the distance between two con- 
secutive collisions, which are as renewal as jumps, but 
are not necessarily jumps. At the collision occurrence the 
system has to decide, through a fair coin tossing prescrip- 
tion, whether to jump to the other state or to remain in 
the original state. This distribution is properly normal- 
ized, and the parameter Tg, making this normalization 
possible, gives information on the lapse of time necessary 
to reach the time asymptotic condition at which i^sij) 
becomes identical to an inverse power law. 

To express the effect of a perturbation on the system, 
we limit ourselves to assuming that 



rP(t) = T s (o)(i±eM*)), 



(3) 



where e is the perturbation strength and £p(i) denotes 
the perturbation signal. A concise account of the dy- 
namic arguments of Refs. [fj Q] is as follows. In the 



2 



limiting case e — * 0, a collision occurring in the "on" 
("off') state at time f', earlier than the collision occur- 
ring at t > t' in the "off" ("on") state, generates the 
response function given by the f'-aged waiting time dis- 
tribution ips(t, t')\ This is the probability density for the 
unperturbed system, prepared at t = 0, to undergo a 
collision at time f, with no collisions happening in the 
interval (f',f), given that observation begins at time t' . 
Thus, we have |E03 



n(t) = <&(*)> = e [ dt'ip s (t,t')Hp(t'). 

Jo 



(4) 



Here H(t) = Pi(t) — p^if), with pi(t) and Pa(^) denoting 
the probability that at time f the system is found either 
in the state |1) or in the state \2). 

Let us explain now how to realize the complex pertur- 
bation £p. The simplest possible way is to adopt as a 
perturbation the signal £p(f) produced by another two- 
state NPR system. The single realization £p(f) is the 
perturbation signal, to which the system S has to re- 
spond. Let us prove first that the perturbed signal £s(f) 
inherits the complexity index \ip. Let us assume that we 
have at our disposal an ensemble of perturbation signals 
£p(f), all of them properly prepared at time t = 0, and 
let us study the double average ((£s(f)))- We realize the 
perturbation preparation by setting all the perturbation 
signals at the beginning of their sojourn in the "on" state. 
It is known 113) that in this case 



<M*)> = = 



Hp— l 



Thanks to the linear nature of Eq. J3J we obtain 



(5) 



(n(t)> = «&(t)» 



dt'^ s (t,t')^ P (t'). (6) 



We establish the time asymptotic properties of (11(f)), 
Eq. |[BJ, by adopting the Fourier expansion of ^(f) = 
J duj^p(uj)exp(—iLut)/(2n), which allows us to use the re- 
sponse to a monochromatic perturbation of Ref. 0, ^| : 
Using this result we obtain that for f — » oo, 



<n(t)> oc 



r( M -i) 



We note that for uj -> £ p (cj) oc w^" 2 . Thus, it 
is straightforward to prove that for f — > oo (11(f)) oc 
l/£' ti>-1 ) with the condition 



Ms + Mp > 3, 



(8) 



whose violation would make the integration over to di- 
vergent. Moving along these lines we explore the ranges 
1.5 < fig < 2 and 3 — pis < Mp < 2, fitting the require- 
ment of Eq. JSJ , and we show 



(H(f)) _ fci(Ms,/jp) ^Cms.Mp) 

e fMP-1 f/iP+l-/iS : 



where 
&i(ms,Mp) 



siii7r^s r(2 - Ms)r(l - + Ms) 

7T r p - MP r(3 - a<p) 



(10) 



For concision's sake we do not write the somewhat ex- 
tended explicit expression of ^(mSiMp); which is not 
used by the heuristic theoretical prescriptions of this 
letter, the only property worth of mention being that 

Mms,Mp) > °> if MS < MP an d k 2 (ns,^p) < 0, if 
Ms > Mp- These analytical and exact results prove that 
the system inherits the complexity index fip of the per- 
turbing signal. 

In the case /ip < /is it is possible to establish a reason- 
able connection between the statistical mean of Eq. , 
expressed in the asymptotic time limit by Eq. @, and 
the response to a single realization of perturbing signal. 
If the system is fast enough as to adapt itself to the per- 
turbation signal £p, as it happens in the case of ordinary 
stochastic resonance 0, 0] , the mean response should re- 
produce the perturbation signal £p. We make the more 
cautious assumption that the mean average adapts itself 
to the signal £p/s(f), which is is very close to £p(f), and 
fits the requirement 

<£p/s(*)> = &>(*)> = *p(*)- (ii) 

Thus, we make the heuristic proposal 

11(f) - e[l - ns{t)]A(iis, Mp)£p/s(*)- (12) 

The function irs(t) denotes the learning function, whose 
time asymptotic property, on the basis of the response to 
constant perturbation [y, [lj] , is assumed to be 



ir s {t) oc f s " 



(13) 



Using Eq. (|llf) . the average over the single realizations 
of Eq. 03) yields 



(11(f)) = e[l - ns(t)]A(fx s , Mp)*p (*), 



(14) 



which is compatible with the theoretical prescription of 
Eq. 10. Thus, using Eq. (|10|) and the asymptotic time 
expression for Eq. JSJ, we conclude that for fis > (J-p the 
amplitude A is given by A(ns, Mp) — ^i(MS) I jl p)/T p if • 
We now have recourse to a numerical treatment for 
the twofold purpose, (i), of establishing the whole time 
evolution of (11(f)) from t = to f = oo, and not only 
the asymptotic time behavior of Eq. 10, and, (ii), of 
establishing the single realizations 11(f). We consider a 
coordinate y moving within the interval / = (0, 1] with 
the equation of motion 



y = a±{t)y* 3 , (15) 

where z s = Ms/(ms - 1) and a±{t) = (ms - l)/Tf{t). 
The motion of y within the interval / corresponds to 

(9) 

the system either in the "on" or "off" state. When the 
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coordinate y reaches the border value y = 1 a new ini- 
tial condition within the interval / is randomly chosen 
with uniform probability. Then we toss a coin to decide 
whether to adopt a + (t) or a„(i) for the next phase of 
this dynamic process. This is the dynamic prescription 
used in 0] to yield the time dependence of Ts, Eq. J3J). 
Note that Eq. © is equivalent to 



n(t) = e / dt'Mt,t)tp(t'), 

Jo 



with 



t + Tp 



/ip— 1 



(16) 



(17) 



Thus, the single trajectories £s(t) are determined by run- 
ning Eq. I|15|) with the time dependent prescription of Eq. 

the proper random back injection and the choice of 
Eq. 1|17|) for the perturbation signal. This allows us to 
check the heuristic prescription of Eqs. (|12|) and iQ . In 
Fig. we show (n(£)) corresponding to /is = 1.6 and 
/ip = 1.5. We see that Eq. (|14fl affords a good fitting 
of this condition. Notice that the adoption of this fit- 
ting prescription allows us to determine A and compare 
it with fci. In Fig. we compare the numerical value of 
A to the theoretical prediction of Eq. 10. We see that 
the agreement between numerical and theoretical result 
is compatible with the statistical accuracy of the numer- 
ical treatment and it becomes worse when /ip > /is, as a 
consequence of /c2 becoming positive and so contributing 
to the slowly decaying tails. 




FIG. 1: Response (TL(t)) for /i P = 1.5 (T P = 1) and /i S = 1.6 
(10 8 systems with Ts = 1, coupling e = 0.1). We assign to 
the learning function the form n s (t) = B 2 " MS /(B + t) 2 ~^ s . 
To fit the data (solid line) we assign to the fitting parameters 
A and B of the overall fitting function 114H (dashed line) the 
values A = 0.80 (in agreement with fci(1.6, 1.5) = 0.798) and 
B — 0.566. Insert: the same in logarithmic scale. 

As to the single realizations, we have to create proper 
perturbation signals, which are complex and do not have 
abrupt jumps that may conflict with the assumptions 




FIG. 2: The amplitude of tA as a function of fip, with Tp = 
1, /is = 1.6, Ts — 1, and e = 0.1 (average over 10 s systems). 
The dots indicate the result of the numerical fitting procedure 
illustrated in the text. For fip < us the procedure rests 
on Eq. 11411 . For fJ,p > Us the fitting function adopted is 
A/t mup ~ l , with t > 100. The solid line is the theoretical 
efci(1.6, /ip). 



made in Ref. 0, Q] to derive the linear response theory 
of Eq. Q}. Thus, we assume that in the natural time 
scale the perturbation is 



£,p{n) — cos(wn). 



(18) 



To make this regular perturbation become complex we 
set first the condition u> <^1. Then, we turn the natural 
time n into the continuous time t with the prescrip- 
tion t(n + 1) = t(n) - 
from the distribution 



r„, where r„ is a number drawn 



V'p(t) = ipp - 1) 



rpHP- 
1 P 



(r + T P y 



(19) 



with £,p(t) = cos(ujt n ) for £„ +1 > t > t n . This procedure 
is inspired to an approach currently used [16i] to derive, 
for instance subdiffusion, by subordinating it, along the 
lines of Continuous Time Random Walk to ordinary dif- 
fusion. It is straightforward to prove that in the time 
asymptotic limit {^p(t)) has the same power law as Eq. 
©. We follow the directions of Refs. 0, El 

and we 

create many perturbed signals £s(i) by means of the pre- 
scription of Eq. . We make the ensemble average on 
a very large number of responses (2 • 10 6 ) and we get the 
results illustrated in Fig. 

The numerical results of Fig. [3] make it possible for us 
to make the heuristic prediction for the response n(t) in 
the case [is < /ip, 



where 



n(i) - e£ s /p(t)R{t), 



(20) 



(21) 



The signal £,s/p(t) is again very similar to the perturba- 
tion £p(i). Also in this case a fluctuating departure of 
£,S/p{t) from £p(i) is admitted. This fluctuation can be 



4 



stronger than the fluctuation of £p/s(i) from £p(i). 
this case we make the request 



In 



(&/p(*)> - (&(*)> - *s(t) 



la 



Ms — 1 



(22) 



Thus, averaging Eq. and using Eqs. and l|2*2*|l . 

we recover the same time asymptotic result as Eq. 
and the constraint (IT(i)} oc l/t tJ/F ^ 1 for f — » oo. 




5000 10000 15000 20000 25000 30000 35000 40000 




5000 10000 15000 20000 25000 30000 35000 40000 
t 

FIG. 3: (A): U(t) for up = 1.8, fis = 1.9. The two lines are 
the numerical response (solid) and, according to the modified 
heuristic prescription (1 1 2t . the signal £p multipied by eA with 
A = 0.34 (dashed). (B): Same for u P = 1.9, fi s = 1.9. Here 
the fitting constant A = 0.37. (C): fip = 1.9, fis = 1.55. 
The two lines are the numerical response (solid) and, accord- 
ing to the modified heuristic prescription 121 H . the signal £p 
multipied by 0.39e(wt) MS " MJ> (dashed), where u = 0.01 is the 
frequency of the natural time harmonic function I18H . 

Fig. prefers to the condition (uj << 1) creating the 
same asymptotic behavior as ^(i) of Eq. JSJ after an 
extended intermediate asymptotic regime, thereby pre- 
venting us from adopting a quantitative comparison with 
Eq. 10 • However, in a qualitative agreement with Eq. 
©; Fig. shows that for fis > fJ-p the response inten- 
sity is weaker than at the matching condition, and that 
for fjbp > /is the response dies out. This proves the CM 
effect, going beyond the theoretical results of Eq. © 
that might generate the wrong impression that the re- 
sponse intensity increases monotonically from fip = 1.5 
up to fip = 2. It is worth pointing out that due to the 
numerical results of Figs. [21 and [3] as well as to the ex- 
act asymptotic time prediction of Eq. for /ip > fis 
we have to replace the heuristic prescription of Eq. 112|) 
with the heuristic prescription of Eq. (|2U[) . which is, in a 
sense, equivalent to the response intensity falling to zero. 
Note that we have modified the heuristic prescriptions 
of Eqs. ljT2|) and (|2*U|l by replacing both £s/p(i) in the 



former case and £p/s(i) in the latter case, with £p(t). 
The agreement between the numerical results and these 
modified heuristic prescriptions is remarkably good. 

In conclusion, this letter shows that a complex system 
with fig < 2 is sensitive to complex perturbations with 
/ip < 2, so that the mean value (H(t)) inherits the per- 
turbation power index fip. The maximum response is re- 
alized through the matching condition fis — fip. If fip < 
(is the response intensity is given by ki((i s , fip)/Tp P ~ 1 
and gets the maximum value when the matching condi- 
tion is realized. When fip > fis the response intensity 
dies out. There are good reasons to believe that with 
A*S > Mp the function £,p/s{t) becomes identical to £p(i) 
thereby yielding a sort of generalization of the well known 
phenomenon of stochastic resonance 0, ^| . 
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